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Large-scale calculations of the El strength are performed within the random phase approximation 
(RPA) based on the relativistic point-coupling mean field approach in order to derive the radiative 
neutron capture cross sections for all nuclei of astrophysical interest. While the coupling to the single- 
particle continuum is taken into account in an explicit and self-consistent way, additional corrections 
like the coupling to complex configurations and the temperature and deformation effects are included 
in a phenomenological way to account for a complete description of the nuclear dynamical problem. 
It is shown that the resulting _Bl-strength function based on the PCF1 force is in close agreement 
with photoabsorption data as well as the available experimental El strength data at low energies. For 
neutron-rich nuclei, as well as light neutron-deficient nuclei, a low-lying so-called pygmy resonance 
is found systematically in the 5-10 MeV region. The corresponding strength can reach 10% of the 
giant dipole strength in the neutron-rich region and about 5% in the neutron-deficient region, and 
is found to be reduced in the vicinity of the shell closures. Finally, the neutron capture reaction 
rates of neutron-rich nuclei is found to be about 2-5 times larger than those predicted on the basis 
of the nonrelativistic RPA calculation and about a factor 50 larger than obtained with traditional 
Lorentzian-type approaches. 

PACS numbers: 24.30.Cz, 24.10.Jv, 21.60.Jz,26.30.Hj 



I. INTRODUCTION 

The investigation of the isovector giant dipole reso- 
nance (GDR) is one of the fundamental problems in nu- 
clear physics and astrophysics. These collective excita- 
tions can be studied experimentally by photoabsorption 
processes and provide essential nuclear inputs for the un- 
derstanding of reaction mechanisms taking place during 
the synthesis of chemical elements in stars. 

Despite new successes of recent large accelerator facil- 
ities, the information we get on the GDR properties is 
still limited to about hundred nuclei in the valley of ft 
stability. All the remaining few thousands exotic proton- 
or neutron-rich nuclei can be accessed by nuclear theory 
only. The aim is therefore to provide the necessary the- 
oretical tools towards a universal description of the nu- 
clear structure phenomena, up to the limit of the nuclear 
stability. 

Over the years, a lot of effort has been devoted to im- 
prove the theoretical description of the main collective 
properties, such as the 7-ray strength function. In par- 
ticular, the random-phase approximation (RPA), based 
on the relativistic mean field (RMF) theory, has proven 
to be a robust theory [J-Q- Pairing correlations for the 
study of open-shell nuclei as well as the self-consistent 
description of the residual interaction Q have also been 
included over the years. 

However, a complete and coherent description of the 
dynamical properties of the nuclei over the entire nu- 
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clear chart remains an open problem and additional cor- 
rections beyond the conventional RPA problem are nec- 
essary, as discussed the present work. 

At first, the coupling to the positive-energy continuum 
has to be taken into account explicitly. This is impor- 
tant, because the alternative approach, namely the dis- 
cretization of the continuum, though simpler, requires 
a basis truncation, which eventually leads to a mixing 
between spurious and physical states of the excitation 
spectrum. In addition, this discretized RPA includes a 
large amount of states which increases considerably the 
numerical effort. 

In fact, the exact treatment of the continuum using 
the Green's function representation not only can solve 
these two problems simultaneously, but also allows for a 
direct calculation of the escape width of the giant reso- 
nance. Although this so-called continuum RPA (CRPA) 
is known to be a powerful approach, it has been rarely 
applied in the past 5 9]. Only recently has a continuum 
QRPA approach been successfully developed within a rel- 
ativistic framework to study giant multipole resonances 
for spherical nuclei [13, [H| • 

Furthermore, many applications, such as nuclear astro- 
physics, require the estimate of the 7-ray strength func- 
tion for deformed nuclei. Recently, microscopic calcula- 
tions with axially deformed RPA approaches |12Hl4j have 
been able to study medium and heavy deformed nuclei 
showing a considerable success in the prediction of the 
GDR properties. In particular they confirm an impor- 
tant effect, namely the splitting of the resonance peak 
into two components. 

Unfortunately, the axially deformed RPA faces several 
difficulties, mostly related to the numerical effort and 
the isolation of the spurious state. For this reason the 
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deformation effect is still treated in all large-scale calcu- 
lations in a phenomenological way. A similar approach 
is followed in the present work. 

The CRPA approach is developed in such a way 
that only one-particle one-hole (lplh) configurations are 
taken into account, so that the coupling to more complex 
configurations, such as 2p2h, is not described in an ex- 
plicit way. These couplings can be taken into account in 
the framework of the relativistic or nonrelativistic time- 
blocking approximation (l5l - [l7l | . However, such a method 
cannot be applied to large-scale calculations yet, essen- 
tially due to the numerical effort it represents, so that 
alternative phenomenological approximations need to be 
considered [3 [l9| . A similar approach is followed in the 
present work, as described in Sec. IIIII 

Regarding calculations of astrophysical interest, it is 
of importance to describe not only the photoabsorption 
processes, but also the reverse radiative particle-capture 
processes. Reaction theory relates the 7-transmission co- 
efficient for excited states to the ground state photoab- 
sorption assuming the giant resonance to be built on each 
excited state. The description of the photo-dcexcitation 
requires however the introduction of a temperature de- 
pendence due to collisions between quasiparticles. As 
traditionally done, a temperature-dependent GDR width 
can be applied and has been shown to improve the pre- 
dictions [2O, nu. 

Finally, as far as pairing correlations are concerned, 
it has been extensively discussed in the past that the 
Hartree-Fock-Bogoliubov (HFB) [22|,[23| or the relativis- 
tic Hartree-Bogoliubov (RHB) approximation [24[ give 
a reliable description of both the ground state and the 
collective excitations of open shell nuclei. A proper treat- 
ment of the pairing correlations is even more important 
in nuclei close to the drip lines where the pairing gap 
becomes of the order of the particle emission threshold. 
However, this approach is not yet available within the 
relativistic CRPA [ll|, so that open shell nuclei are still 
treated within the BCS model applied in the mean field 
as well as the dynamical system. Until an RHB plus 
CRPA is developed, the BCS model can be regarded as 
a useful approximation of the pairing effect on the 7-ray 
strength function. 

The purpose of the present paper is to determine sys- 
tematically the CRPA El strength function for the entire 
nuclear chart within the relativistic mean field approach. 
For practical applications, all the above mentioned cor- 
rections (correct description of the continuum, spreading 
width, deformation effects, temperature effects) are in- 
cluded for a proper prediction of the reaction cross sec- 
tion. The paper is organized as follows. In Sec. [XT] we 
briefly describe the relativistic quasiparticle RPA, where 
the continuum is explicitly included. In Sec. IIIII the 
damping of the collective motion, the impact of defor- 
mation on the GDR and the temperature dependence 
are included in the model and the new strength is com- 
pared with experimental data. In Sec. IIVI a large-scale 
calculation for about 8000 nuclei is shown to give inter- 



esting properties of the nuclear collective motion. The 
proton pygmy resonance of neutron deficient nuclei is ex- 
tensively studied in Sec. |V] Finally, in Sec. |VI] we use 
the final strength to study the radiative neutron capture 
rates for all the 8000 nuclei and report our conclusions 
in Sec. EED 



II. RELATIVISTIC CONTINUUM QRPA 

As in all the relativistic models, the nucleons are de- 
scribed as point like Dirac particles, which move inside a 
nuclear mean field. This mean field is phenomenological 
and can be described by the exch ang e of effective mesons, 
as is done in the Walecka model [251 ] or by using the zero 
range or point coupling terms. This latter description of 
the RMF introduced in the early 1990s [H H3 has since 
been improved substantially and is used in our work to 
describe the ground state and collective properties. This 
approach has the main advantage of simplifying consid- 
erably the numerical calculations without loosing in pre- 
dictive power as compared to the meson exchange mean 
field. 

In this work, we use the point coupling Lagrangian 
PCF1 [28| expanded in powers of the nucleon scalar (S), 
vector (V) and isovector-vector (TV) fields: 



C = 'ip(i"f ■ d — m)rp 



7 (i-^) , 
eip-f ■ ip, 



(c=S,V,TV) (1) 



where the Dirac vertices r c have the explicit form 
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(2) 



For the PCF1 force, the index c takes just seven val- 
ues, namely one scalar S, three isoscalars V and three 
isovector- vectors TV. The gradient terms S c (d^T c ^fj) in 
the Lagrangian C are used to simulate the finite range 
of the nuclear interaction and consequently improve the 
description of the nuclear surface properties. The energy 
functional deduced from Eq. (fTJ) is sufficient to describe 
the nuclear ground state, the single-particle Hamiltonian 
being given by: 



Spa/3 



(3) 



To describe the nuclear response to an external field F 
induced by a particle or by photoabsorption, the strength 
function 



S(E) = --lm V F* p R aPa ,p,{E)F a , fi , (4) 

af3a'P' 
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needs to be estimated. The Greek indices a, (3 indicate 
the various degrees of freedom of a nucleon (r,s,t,d), 
where s is the spin, t the isospin coordinate, and d = 1, 2 
labels the large and small components. If we consider the 
case of a weak external field and thus a small amplitude 
variation of the static solution, the response function 
Rafia' 0' {E) can be deduced from the linearized Bethe- 
Salpeter equation 



Ra/3a'l3'(E) 



E 



(5) 



The residual interaction V? s , s , of Eq. (0 is connected to 
the static problem via the second derivative of the energy 
functional 
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P h 
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This expression of , s , includes in addition to the lin- 
ear time-like fields, the space-like terms, the Coulomb 
part and the rearrangement terms, as a result of the sec- 
ond derivation 11| . In this way, a self-consistent solution 



is ensured, which is of particular importance for the de- 
scription of exotic nuclei for which no data exists. 

In coordinate representation, the indices a,/3, . . . in 
Eq. ([5]) are abbreviations for the "coordinates" 1 = 
(ri, S\, t\, d\). The simplicity of the point coupling La- 
grangian allows us to express the effective interaction of 
Eq. ([6]) as a sum of separable terms, i.e 



V* h (l,2) =J2Q*c {1) (r)v c (r)QW(r) 



(7) 



with the local channel operators Qc (r) defined by 

g(i) (r) = ( _ )So fc^)fWF i (0 1 )- (8) 
rr\ 

These single-particle operators are characterized by the 
channel index defined in the spin-isospin space (c) as well 
as the coordinate space (r). Therefore, assuming a coor- 
dinate mesh r of 50 points and having c = 7, the size of 
the interaction matrix is typically smaller than 350 x 350 
which significantly reduces the numerical effort . 

Finally, by inserting the effective interaction of Eq. ([7]) 
into the Bethe-Salpeter equation (0, we get the expres- 
sion 

K cd (E) = K° c , (E) + Kc>> {E)v c „ {r)Tl c n cl (E), (9) 



where the reduced response function is given by 



n cc ,(E) = q* c ^r(e)q [ J 



l(2) 



(10) 



Equation ^ has the same formal solution as Eq. ([5]), but 
is simpler to solve. At the end, one has to calculate the 



free response function 71® C , (E) which corresponds to the 
basic quantity of the CRPA approach. 

According to 0, EH, in order to treat the coupling 
to the continuum exactly, the nonspectral representation 
T^cont needs to be estimated from 



Kont = J2 v k( k (r)\Qc9(r,r';E-E k + X) (11) 



+g(r',r:-E-E k + X)Ql,\k(r')) 



The Green's function is defined here as 



g(r 1 ,v 2 ;e K ) = — j 



Ly^/ KM)«(r')l r>r' 

K(r))«(r')| r<r>, ^ 



where the Wronskian W acts as a renormalization factor. 
The two-dimensional spinors \u K (r)) and \w K (r)) are the 
regular and irregular wavefunctions of an excited state k. 

The nonspectral method [Eq. (TT2|) ] has several advan- 
tages. First, for energies above the e K = continuum 
limit, the function \w K (r)) is a plane wave with complex 
values. In other words, for E > A (which accounts for all 
excitation energies above the particle emission threshold) 
the strength function [Eq. |4])] is always nonzero. This 
finite distribution of S(E) is a direct source of informa- 
tion for the escape width of the giant resonance [l(| ■ 
Furthermore, all levels embedded in the continuum and 
included in the Green's function [Eq. (jT2"j) ] arc character- 
ized by the quantum number k and not by the principal 
quantum number n, as done in the discrete approxima- 
tion. In other words, there is no distinction between e.g. 
states 2j>3/ 2 and 3p 3 / 2 - In this way, the sum over a large 
number of unbound states is reduced to a sum of only 
a few k states. In addition, the excitations to the Dirac 
sea (which includes the antiparticle states) are taken into 
account automatically and thus their calculation is re- 
dundant. Both these characteristics can lead not only to 
a substantial reduction of the numerical effort but also 
to the absence of a corresponding cut-off parameter, as 
required in the discrete RPA approach Q. 

As discussed before, for open-shell nuclei, the BCS 
approximation is used at both the RMF and the RPA 
levels in order to take the pairing correlations into ac- 
count. As a consequence, the bound states in the pairing 
active space are described by an occupation factor, i.c, 
as quasiparticles. However, since the continuum states 
above E = remain pure particle states, the bound 
and unbound states are not treated on an equal foot- 
ing anymore. For this reason, the response function of 
Eq. (jlip cannot be properly described within the quasi- 
particlc CRPA model which needs to be extended to ac- 
count for the two-quasiparticle excitation in the pairing 
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active space. This is done by adding the term 

KpW = E T^-{k\\Qo\W)r{k\\Q c '\\k% 
k<k , l + Ofcfc' 



1 



1 
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<E P 

+ 52 i , X (fcllQcllfcOr^HQc'llfcV 

f-ri, 1 + Okk' 



k<k' 



e — sit k' e -\- n 



E — ilk', k E + ilk' 



(13) 



where Ekk' = Ek + Ek> is the two quasiparticlc energy, 
E = E + iT/2 and tt kk ' = E k - e k > - A. In Eq. (13]), the 
first part corresponds to two-quasiparticle excitations, 
while the second part is a correction term to prevent the 
excitations in the pairing active space from being counted 
twice (sec Rcfs. [H [ll[ for more details). Equation (fT3"|) 
does not imply any significant additional numerical effort 
since the sum applies only in the pairing active space of a 
few quasiparticlc excitations. Consequently, a quasipar- 
ticle CRPA requires that the extended response function 
is properly treated only after the use of 



TZ°(E)=TZ° cont (E) +K° 2 (E) 



(14) 



in Eq. ©. The cross section for the photoabsorption 
processes is then derived from the strength function of 
Eq. (Q| according to 



1 PrtrSpZ 

&(E) = -^-^ES RPA (E) [fm 2 ] 
9fi c 

= 4.022 ES RPA (E) [mb]. 



(15) 



III. CORRECTION TO THE CRPA STRENGTH 

A. Damping of the El strength 

Although the escape width is estimated microscop- 
ically within the CRPA, its value remains small, espe- 
cially for heavy nuclei. This is mainly due to the high 
centrifugal and Coulomb barriers, which prevent the ex- 
cited nucleon from escaping. For that same reason, in 
light nuclei, where the Coulomb and centrifugal barri- 
ers are smaller, the escape width dominates the total 
width 0. 

The remaining part of the total width corresponds to 
the coupling of the resonance to more complex configu- 
rations (2p-2h, 3p-3h, etc.) and its contribution to heavy 
nuclei is large. One way to deal with the damping width 
is to calculate it in the second-RPA (SRPA) framework 
[l8j or by particle vibration coupling [l5l - [l7| . 



Formally, self-energy insertions on particle and hole 
states spread the resonances and shift their centroids. In 
a first approximation, this can be done by folding the 
RPA strength with a Lorentzian function 



f (E j?) = I T (E)E* 

' k[E 2 -{E>- A(E') 2 } 2 + T(E) 2 E 2 

so that 



, (16) 



S(E)= dE'S R PA(E')f L (E,E'). (17) 

JQ 

Here E' is the excitation energy of the CRPA response, 
and A(E) and T(E) the real and complex part of the self- 
energy, respectively (l8| . The energy- dependent width 
T(E) can be calculated from the measured decay width 
of particle (j p ) and hole (7/J states [l8| 

r(E) = ^j\k[ lp (E) + lh (-E)](l +45). (18) 

This empirical way of determining T(E) has the advan- 
tage of including, in principle, contributions from the 
excitation beyond 2p-2h [ljj. The resulting resonance 
width can be compared with experimental data, such 
as photoabsorption cross sections. Finally, the real part 
A(E) of the self-energy is obtained from T(E) by a dis- 
persion relation [l~8[ , where the interference factor Coy is 

allowed to differ from C sl i affecting the Lorentzian width 
of Eq. p|l. 

This method, though empirical, avoids important 
SRPA difficulties [29| related to the large deviation of the 
GDR energy and the need to renormalize the residual in- 
teraction in order to deal with the spurious components 



B. Comparison with the experimental GDR 

While the GDR of spherical nuclei can be well de- 
scribed by a Lorentzian-type function (Eq. ITB|) , the situ- 
ation for nuclei with a deformed ground state is different. 
In particular, in the case of deformed nuclei with axial 
symmetry, the GDR is known to split into two major 
resonances which correspond to oscillations of protons 
against neutrons along the principal axes, with frequen- 
cies inversely proportional to the length of these axes 
[l3l |. In the phenomenological approach, the GDR peak 
energy is split according to the following rule [33| 



E 



(i) 

GDR 



(19) 



4dr /4dr = 0.9117? 



0.089, 



where rj is the ratio of the diameter along the axis of 
symmetry to the diameter along an axis perpendicular 
to it. Consequently, the folding process of Eq. (fTTf takes 
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place twice for each energy, with the two Lorentzians cor- 
responding to the two energy peaks of Eq. ([T9"]) , leading 
to the final strength 



S(E) 



dE'S RPA (E') 



1 



f£>(E,E') + ?-tf>(E,E> 



(2) 



(20) 

This procedure is applied to all nuclei, the parameter rj 
being derived from the HFB mass model of Ref. [34[ ■ 

The two interference factors C ST and C^ T related 
to the energy shift and the resonance width of the El 
strength, respectively, are calculated in a way similar to 
the parameter set of the Lagrangian, i.e. they are ad- 
justed through a minimization procedure to reproduce 
at best the measured GDR energy and width of experi- 
mentally known nuclei. In particular, the rms deviation 
factors 



f(E) 
i rms 



i 



N, 



£ [E z G B R (th) - E* GDR (exp)] 



1/2 



(21) 



T) . 

ms 



l N - 

[ r GDR.(^) - r GDR.(e^)] 

G 2 — 1 



1/2 



(22) 



are used for the GDR energies and the GDR widths, 
respectively. The best fit to experimental data is ob- 



tained for C< 



(E) 
ST 



-0.85 and C, 



f r ( il = 0.66 MeV and /£? B 



ST 



-0.59, leading to 



1.81 MeV. This result is 
illustrated in Fig. [TJ where the GDR energy and width 
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FIG. 2: (Color online) Comparison between the experimen- 
tal photoabsorption cross section approximated by a simple 
Lorentzian curve (black dashed line) and the CRPA predic- 
tions obtained with the PCF1 force (red solid line), for nine 
representative spherical nuclei given by (Z,A). All cross sec- 
tions are normalized to a peak cross section of unity. Experi- 
mental data are taken from [3~D . I32I ] 



are given for about 80 nuclei and compared with existing 
experimental measurements. 

In Fig. [2] we illustrate the El photoabsorption cross 
sections (normalized to unity) for several spherical nu- 
clei all over the nuclear chart. We see that globally, the 
GDR widths and energies are in good agreement with ex- 
perimental data (represented by a Lorentzian function) 
[3ll l32j ]. When the deformation effect becomes impor- 
tant, it can be seen in Fig. |3] that the folded CRPA 
strength [corrected with the prescription of Eq. (TIT))) ] pre- 
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FIG. 1: (Color online) (a) Comparison between the experi- 
mental GDR energies for spherical and deformed nuclei (black 
squares) and the GDR peak energies (red cycles) obtained 
after folding the CRPA distribution and including the defor- 
mation effects, (b) Similar comparison between experimental 
and calculated GDR widths. 



FIG. 3: (Color online) Comparison of the CRPA predictions 
(solid red line) with the experimental photoneutron cross sec- 
tions (black dashed line) for several nuclei with well-deformed 
ground state. The experimental data are taken from [35M401] 
for the nuclei 150 Nd, 168 Er, 172 Yb, 181 Ta, 237 Np and 238 U 
respectively. 
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diets rather well not only the GDR width, but also the 
splitting into the two observed peaks. 



C. Temperature dependence 

As already discussed, when dealing with practical ap- 
plications and more precisely radiative captures, the de- 
excitation strength function needs to be estimated. For 
low-energy incident particles, a reliable description of the 
tail of the GDR close to the particle threshold is required. 

The excitation of a nuclear state in this energy range 
is followed by a subsequent decay to the ground state 
which takes place via multiple intermediate states of fi- 
nite lifetime. The observation that in many nuclei the 
deexcitation 7-strength at E 



— ► is nonzero [2l| sug- 
gests that the level density in this region is high and that 
the El str eng th function depends on the nuclear temper- 
ature (43l . |46| . To provide a qualitative agreement with 
7-decay measurements, a temperature dependence is tra- 
ditionallyintroduced in the expression of the GDR width 
SI El El, i.e., 



T'(E,T) = T(E) 



1 



rp2 
GDR 



E l 



a47T 2 T 2 E GDR 
E + 5 



(23) 



where T is interpreted as the nuclear temperature of the 
final state and is estimated by Tf oc WEf. The con- 
stant parameter 5 — 0.1 MeV is introduced to ensure a 
nondivergent width at E — > (2lj . 

Available data on thermal capture measurements from 
Mg to U have been used as a benchmark for the validity 
of the model and also as a way to derive the parameter a. 
In Fig. 0] the corresponding strength functions fj(E) are 
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FIG. 5: (Color online) CRPA 7 -ray strength of 96 Mo, 98 M 
and 144 Nd, obtained with a T- independent (blue dotted line) 
and T-dependent (red solid line) width. The results are com- 
pared with experimental data for the Mo [ill . |4^ | and 144 Nd 
isotopes mEEl- The CRPA calculations corresponds to a 
temperature T = 0.50 MeV. 



shown at energies ranging between 4 and 8 MeV depend- 
ing on the corresponding binding energies. It is found 
that the value a — 2.30 gives the best overall agreement 
with the experimental strength. We see that the temper- 
ature dependence improves the prediction, as compared 
to the T ~ strength. 

The importance of the temperature dependence can be 
even more apparent for nuclei with enhanced El strength 
at low excitation energies, as for instance in the deexci- 
tation strength of 96 Mo, 98 Mo or 144 Nd (Fig. |5J. After a 
proper rcnormalization of the position of the GDR, the 
experimental strength at E ~ is found to be well re- 
produced by the convoluted strength function with the 
temperature-dependent width of Eq. 02 
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FIG. 4: (Color online) Comparison of the T-dependent and T- 
independent low-energy .El-strength functions deduced from 
the CRPA with the experimental compilation [32] including 
resolved-resonance and thermal-captures measurements, as 
well as photonuclear data for nuclei from 25 Mg up to 239 U 
at energies ranging from 4 to 8 MeV. 



IV. LARGE SCALE CALCULATIONS 

Large-scale calculations based on the extended CRPA 
have been performed for all 8 < Z < 110 nuclei lying 
between the proton and neutron drip lines, i.e. some 
8000 nuclei. For all the even-even nuclei, the CRPA 
El distribution is folded by a Lorentzian-type function 
which properly takes into account the effect of the in- 
terference factors, the deformation and the temperature 
dependence, as described in Sec. IIIII Interpolation tech- 
niques are used to calculate the strength of the odd-^4 and 
odd-odd nuclei. This interpolation is based on a simple 
linear averaging between the two closest even-even nuclei, 
with respect to the GDR energies, i.e for an odd nucleus 
(Z, N), the El strength is given by 



S^ N \E) 



S( N -V(E + A/2) + S( N +V(E - A/2) 



where the quantity A = E, 



(JV-l) 
GDR 



-E, 



(JV+1) 
GDR 



corresponds to 



the difference in the GDR ccntroid energies between the 
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FIG. 6: (Color online) CRPA centroid energies (a) and widths 
(b) of even-even A nuclei with 8 < Z < 110 and 10 < N < 
240. 



nuclei {Z,N — 1) and (Z,N + 1). Despite the empirical 
nature of this method, experimental results on low-lying 
El strength j4?| and GDR data [32[ suggest that it can be 
adequately used in medium and heavy nuclei. A proper 
microscopic treatment of odd nuclei has been developed 
[48l ] , but remains numerically too heavy to be applied to 
large scale calculations. 

In Fig. [6l the GDR centroid energies and widths are 
plotted with respect to the neutron number N for all 
even-even nuclei. It can be seen that the GDR en- 
ergy decreases with the atomic mass while the GDR 
width reaches local minima in the vicinity of magic num- 
bers, thus being sensitive to shell effects. For extremely 
neutron-rich nuclei with N > 200, the GDR width ap- 




pears to remain rather constant around 8 MeV. 

Apart from the GDR properties which dominates the 
photoabsorption cross section, a second weaker mode ap- 
pears at low energy, especially for nuclei with a large neu- 
tron excess. This pygmy dipole resonance (PDR), which 
has been extensively studied in the last decade (see, e.g., 
[49l - l5ll |). is found at low energies, typically around 9 MeV 
for nuclei around the valley of ft stability and at even 
lower energies for nuclei close to the drip lines (Fig. [7]). 
Since this excitation mode lies in the low-energy tail of 
the GDR, it is expected to have a significant impact on 
the determination of the radiative neutron capture rates 
of exotic neutron- rich nuclei [49|, HO, HH . 

In the upper panel of Fig. [5] the energy of the pygmy 
resonance is shown for all even-even nuclei up to N = 
240. We see that the PDR energy decreases as we move 
to more neutron rich nuclei, a pattern similar to the 
one characterizing the GDR. This decrease is however 
slower and for superheavy nuclei, -Epdr remains rather 
constant. In Fig. IHb), the integrated PDR strength 
relative to the full El strength is given as a function 
of the neutron number N . The PDR strength is ob- 
tained by integrating over the energy region around the 
PDR peak energy and up to about 10 MeV, where the 
GDR contribution dominates. The behavior of the PDR 
strength shown in Fig. [5] demonstrates the shell effects 
affecting the pygmy mode; close to the neutron magic 
numbers N =20, 50, 82 and 126 a significant decrease of 
the PDR strength is observed, while far from the magic 
numbers, this strength can become as large as 15% of the 
total strength. Finally, for the neutron-rich nuclei with 
N > 200, the PDR strength remains constant at about 
5% of the GDR strength. 
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FIG. 7: (Color online) CRPA El strength function for the Sn 
isotopes with A = 116 to A — 160 by steps of 4 as a function 
of the energy. 



FIG. 8: (Color online) Same as Fig. [6] for the position of the 
PDR (a) and the ratio of PDR strength to the full El strength 
(b). 
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V. PROTON PYGMY MODES 

The nature of the PDR in neutron-rich nuclei, as dis- 
cussed in the previous section, refers to the well-accepted 
picture of a neutron skin vibration with respect to a 
neutron-proton core. 

In principle such a picture can be generalized to 
neutron-deficient nuclei for which an oscillation of the 
weakly bound proton skin could take place against the 
isospin-saturated proton-neutron core. This idea is of 
course not new and has been studied in the past, though 
essentially on a theoretical ground [H, HH . 

While many neutron-rich system exist, neutron- 
deficient nuclei are restricted to Z < 50 elements, 
the Coulomb barrier preventing the creation of heavier 
proton-rich systems. In light neutron-deficient nuclei, 
the proton excess remains rather small, so that a proton 
PDR is expected to be rather weak, at least in compar- 
ison with the neutron PDR on the neutron-rich side of 
the valley of /3 stability. In addition, the El strength in 
light nuclei is quite spread, blending a possible proton 
PDR contribution. 

On the basis of our theoretical approach, we can iden- 
tify and study the properties of the proton PDR. As 
shown in Fig. [9l the neutron-deficient 28 Ar and 34 Ti iso- 
topes reveal an enhanced El strength around 6.9 MeV 
and 7.4 MeV, respectively, which could be identified as 
a proton PDR. The first indication for this is that the 
low-lying strength disappears, as we reach the N = Z 
isotopes, i.e. for 36 Ar and 44 Ti. Furthermore, the proton 
and neutron transition densities (lower panels of Fig. [S]) 
clearly show no contribution of any neutrons in the sur- 
face region; instead, we observe the existence of a proton 
skin which oscillates against the proton-neutron core, as 
suspected in the first place. 
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N 

FIG. 10: (Color online) Ratio of the PDR to GDR strength 
with respect to the neutron number N for the even-even 
proton-rich nuclei, from the proton drip line up to the N = Z 
isotopes 

In Fig. [TU1 we estimate the proton PDR strength rel- 
ative to the GDR strength for the even-even neutron- 
deficient nuclei with Z < 50. The quantity ttipdr cor- 
responds to the integrated El strength in the 5-10 MeV 
region. For each isotopic chain, we observe a rapid de- 
crease of the PDR strength, as the mass number increases 
towards the N = Z line or the magic numbers N = 20, 
28 or 50. This is a common behavior for all light and 
medium nuclei, except for some specific chains around 
Z ~ 32. This deviation from the general trend is also 
found for the neutron PDR strength for nuclei around 
N = 50, N = 82 and N ~ 144, as seen in Fig El 



VI. RADIATIVE NEUTRON CAPTURE RATES 



5 0. 
„ § 0. 

go. 



1 1 1 1 1 

/ 


1 1 


- — N=10 




- --■ N=14 i 




_ — ■ N=18 j^, 




" / 






Z=18 " 







5 10 15 20 
E [MeV] 



Ar- 



2 4 



r [fm] 



3- 
2- 
l- 
0- 
-1 L 
-2- 
-3- 

_i_LJd -4L 
10 12 



6.9 MeV 

I I L 



Z=22 | 



10 15 20 
E [MeV] 




Ti 



7.4 MeV 



_l 1 L 

4 6 
r[fm] 



8 10 



FIG. 9: (Color online) (Upper panels) El strength of the light 
Ar and Ti isotopes, showing a pronounced proton skin, as 
compared to the N = Z isotopes (dot-dashed lines). (Lower 
panels) 28 Ar (left) and 34 Ti (right) proton and neutron tran- 
sition densities at the PDR. 



The Maxwellian-averaged radiative neutron capture 
rates of astrophysical interest (sec c.g the review [HH 
on the r-process nucleosynthesis) are estimated within 
the statistical model of Hauser-Feshbach at a typical 
stellar temperature of 10 9 K, making use of the TALYS 
code (5(| [E3|- This version benefits in particular from 
an improved description of the nuclear ground state 
properties derived from the nonrelativistic Hartrce-Fock- 
Bogoliubov (HFB) method 58| . as well as from a nu- 



clear level density prescription based on the combinato- 
rial model [sH]. The direct capture contribution as well 
as the possible overestimate of the statistical predictions 
for resonance-deficient nuclei are effects that could have 
a significant impact on the radiative neutron captures 
by exotic nuclei [49| . but arc not included in the present 
analysis. 

In Fig. 1111 we compare the radiative neutron cap- 
ture rates (av) for the Sn isotopes obtained with 
our T-dependent CRPA with both the nonrelativistic 
HFB+QRPA predictions [2(| and the widely used phe- 
nomenological generalized Lorentzian model (GLO) of 



Ref. [43| . Although all three models lead to reaction rates 
within a factor of 2 for Sn nuclei below TV = 82, this is 
not the case for exotic neutron-rich nuclei, for which the 
CRPA model gives rise to radiative neutron capture rates 
about 10 to 60 times larger than those obtained with the 
GLO model. This can be explained by the presence of a 
strong low-lying CRPA strength for neutron-rich nuclei, 
as shown in Figs. [7] and [5] 

A similar, though less pronounced, pattern is observed 
with respect to the nonrelativistic HFB+QRPA approach 
[20j . In fact, in the region between 132 Sn and 142 Sn, re- 
action rates larger by a a ratio of about 5 is found with 
our CRPA approach. This is due to the fact that the 
neutron threshold for these nuclei lies at energies close 
to the pygmy resonance, which is generally more pro- 
nounced in the relativistic model. In contrast, heavier 
isotopes, i.e. 146 Sn or above have a neutron threshold 
below 3 MeV and are less affected by the PDR, which, 
although stronger, lies at higher energies. For this reason, 
both models tend to give rise to similar reaction rates for 
the most exotic neutron-rich Sn isotopes. Calculations 
based on a relativistic [6(J or nonrelativistic [l7| QTBA 
description of the Sn isotopes show a similar behavior 
when compared to the HFB+QRPA approach [2(|, al- 
though the low-lying strength obtained within the time- 
blocking approximation appears to be more fragmented 

The ratio of the CRPA neutron capture rates to 
those calculated with the nonrelativistic HFB+QRPA 
are shown in Fig. [12] for all nuclei with 8 < Z < 92. 
Globally, the CRPA calculation tends to predict larger 
El strength at low energies mainly due to the lower ener- 
gies at which the PDR dominates. In some very neutron- 
rich nuclei the rates calculated with the CRPA strength 
can be about 10 times larger. Finally note that the In- 
dependence introduced on the CRPA strength in Eq. (|2"3")l 
plays a minor role on the neutron capture rates of exotic 
neutron-rich nuclei since their low neutron separation en- 
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FIG. 11: (Color online) Ratio of the Maxwellian-averaged 
(n,y) rates (at a temperature of 10 9 K) for the Sn isotopes 
obtained with the CRPA El strength to the one using the 
HFB+QRPA [U (squares) or the GLO model g| (circles). 




FIG. 12: (Color online) Ratio of the Maxwellian-averaged 
(n,y) rate (at a temperature of 10 9 K) obtained with the 
CRPA El strength to the one using the HFB+QRPA [U 
for all nuclei with 8 < Z < 92 lying between the proton and 
neutron drip lines. 



crgy implies a low temperature of the dcexciting levels. 



VII. CONCLUSIONS 

Starting from a point coupling Lagrangian with the 
parameter set PCF1 [28[, we have used the relativistic 
continuum QRPA approach to study the El collective 
excitation spectra of all nuclei between the neutron and 
proton drip lines. The RMF equations are solved in the 
coordinate space sclf-consistcntly. For open-shell nuclei, 
the BCS model is applied to treat the pairing correla- 
tions. 

The residual particle-hole interaction used in the RPA 
calculations is derived from the same Lagrangian in a 
fully self-consistent way. A set of two interference factors 



Cc^j} and Cg T ; is additionally introduced in order to ad- 
just systematically the damping width and the centroid 
energy of the GDR energy on all photoabsorption data. 
In addition, the deformation and temperature effects are 
included to account for a complete description of the nu- 
clear dynamical problem. This extended CRPA model 
has been tested for nuclei for which experimental infor- 
mation is available and a satisfactory quantitative agree- 
ment is found. In particular, the predicted GDR peak en- 
ergies globally reproduce experimental data within only 
a few hundred kcVs. 

On this basis, the extended CRPA model is used to 
perform a large-scale calculation of the El strength for 
all nuclei with 8 < Z < 110 lying between the proton and 
neutron drip lines. As far as the collective dipole excita- 
tion is concerned, a clear pattern around magic numbers 
is observed on the GDR width and the PDR strength, 
both being affected by the shell effects. The tempera- 
ture dependence of the width also appears to affect the 
strength in the low-energy region. 

For neutron-rich as well as light neutron-deficient nu- 



10 



clei, a low-lying PDR strength is found systematically in 
the 5-10 MeV region. The corresponding strength can 
reach 10% of the GDR strength in the neutron-rich re- 
gion and about 5% in the neutron-deficient region, and 
is found to be significantly reduced in the vicinity of 
the shell closures. Finally, the neutron capture reaction 
rates of neutron-rich nuclei determined with the CRPA 
strength can be about 2-5 times larger than those pre- 
dicted on the basis of the nonrelativistic HFB+QRPA 
strength and up to about 50 times larger than classically 
determined with the GLO model. 

We have shown that the present approximation can 
described rather successfully the collective properties of 
stable nuclei for which experimental data are available. 
However, for more exotic nuclei, in particular close to the 
drip lines, some of the phcnomenological corrections used 
in the present approach need to be replaced by sounder 



models. More specifically, an improved treatment of the 
pairing correlations can be achieved within the Relativis- 
tic Hartree-Bogoliubov model [HI and the effects linked 
to the spreading of the GDR strength beyond the lp-lh 
excitations as well as the ground state deformation and 
odd number of particles need to be tackled following the 
recent works of [lj HI EE El- 
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